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Abstract 

We compute, within the Schrodinger functional scheme, a renormaliza- 
tion group invariant renormalization constant for the first moment of the 
non-singlet parton distribution function. The matching of the results of our 
non-perturbative calculation with the ones from hadronic matrix elements 
allows us to obtain eventually a renormalization group invariant average 
momentum of non-singlet parton densities, which can be translated into a 
preferred scheme at a specific scale. 
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Physical quantities that need renormalization, such as the coupling constant, the 
quark mass or the matrix elements of operators appearing in the Wilson opera- 
tor product expansion with a non-zero anomalous dimension, are "running" with 
the renormalization scale. The choice of the scale is in general motivated by the 
kinematics of the Green functions involving such renormalized quantities, but 
the final physical predictions of the theory without perturbative approximations 
are independent of such a choice. This leads to the well-known renormalization 
group equations that put the independence on a formal basis. The redundancy in 
a parametrisation of the theory in terms of renormalized quantities and the rel- 
ative renormalization scale can be avoided by considering renormalization group 
invariant quantities, such as the A parameter of QCD or the renormalization 
group invariant quark masses (RGIM). The advantage of the latter choice in non- 
perturbative lattice determinations of the quark mass has recently been stressed 
by the authors of ref. |1| where an essential part of the RGIM programme was 
carried out. 

In particular, the definition of the RGIM, which corresponds, roughly speaking, 
to a running mass at infinite renormalization scale, is free of the renormaliza- 
tion scheme dependence that usually affects quantities renormalized (in a given 
scheme) at a fixed scale. It can hence be evolved back to an arbitrary finite scale 
in a preferred scheme. 

The purpose of this letter is to present a similar calculation for the operator 
that corresponds to the average momentum of non-singlet parton densities. A 
lattice — perturbative and non-perturbative — study of the evolution of such an 
operator has been discussed in refs. M and M to which we address the reader 
for more details about the calculation that we here only shortly summarize as 
follows. We calculate the renormalization constant of the twist-two non-singlet 
operator for the first moment of the quark parton distribution defined by: 

Of = {\T~ l ^)l{, D v} yV(^) - trace terms , (1) 

where {. . . } means symmetrization of the indices. We remark that the technique 
discussed here can be extended to higher moments in an anloguous way. The 
basic ingredient for the reconstruction of the non-perturbative scale dependence 
of the renormalization constants of the above operator is the finite-size step scaling 
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function az defined by: 

Z(sL) = a z (s,g 2 (L))Z(L) , (2) 

where L is the physical length that plays the role of the renormalization scale, s 
parametrizes the step size of the change in the scale, and Z is the renormalization 
constant of the operator, which is defined by: 

R (fi) = Z(l//i)- 1 C barc (l/L) . (3) 

Z is obtained from the Schrodinger Functional (SF) matrix element, (. . . )gp, of 
the operator on a finite volume L 3 T, normalized to its tree levelQ: 

(O barc (l/L)) SF = Z(L)(O trcc ) SF . (4) 

The renormalized operator then satisfies {O r (h = 1/L))sf — (C tree )sF- The 
framework of the Schrodinger Functional [f|, which describes the quantum 
time evolution between two fixed classical gauge and fermion configurations, de- 
fined at times t = and t = T, has been used extensively in the recent literature 
P, |7|, [1[ to calculate non-perturbative renormalization constants of local oper- 
ators. Among the advantages of the method, we only quote the possibility of 
performing the computations at zero physical quark mass and of using non-local 
gauge-invariant sources for the fermions without need of a gauge-fixing procedure. 
In our particular case, we exploit both features. Our observable is defined by D: 



f 2 (xo = L/A) /ff 2 (x = L/4) 



(5) 



tree 



with / 2 given by 

/,(*„) = -a 6 ^e^- z )(^(x)7{i D 2} ^(x)C(y)r^T 3 C(z)> (6) 



and fi by 



r 3 _ ^, ,r 3 



h = -a 12 E (C'(v)-C'(w)C(x)-C(y)) , (7) 



y,z,v,w 



1 In the following we choose T = L. 
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where ( = 5/5ip c and ( = —6/6ip c are the derivatives with respect to the two- 
component classical fermion fields (z/> c and ip c , respectively) at the boundary xq = 
0, while (' and (' are the corresponding derivatives at the boundary xq = T. The 
projection on the classical components is achieved by the projector P± defined by 
|(1 ±70). On the boundaries, the theory possesses only a global gauge invariance 
that is preserved by the quantities defined above. The values of Xq (set to T/4) 
and of the non-zero component of the momentum p x (set to 27r/ L) are both scaled 
in units of L, which therefore remains the only scale besides the lattice spacing a. 
The quantity f\ serves as a normalization factor that removes the wave function 
renormalization constant of the ( fields in order to isolate the running associated 
with the operator in eq. (|1|) only 

The determination of the step scaling function in the continuum has been shown 
to be universal with respect to the lattice action used in ref. ||. From a fit to 
its dependence upon the running coupling constant g 2 , renormalized in the SF 
scheme, we can extract the following "running" step scaling function: 



<7(////i , f(»o)) = Z(l/fx)/Z(l/fx ) (8) 

i.e. the renormalization constant normalized to the one at a reference scale Hq. 
The running operator matrix element at the scale fi, which we denote generically 
by the symbol O, can be defined in terms of the one at scale /io simply by: 

0»(ri=0™Wa(^„ l5 2 W). (9) 

The scale dependence of the renormalized operator just reflects the one of its 
renormalization constant governed by the equation: 

Z(l//i)-7oG?V)), (10) 



dZ(l/fi) , ,, 



dlog(/i) 
from which follows: 

dO ren (/i) 



O ren (/i)-7oG? 2 M). (11) 



d log(/x) 

Following ref. [|IJ but using a slightly different normalization in taking out the 
factor of 2&o, we define, for operators entering the Wilson operator product ex- 
pansion, a renormalization group invariant matrix element: 

OSv = O re » ■ (sT 70/26 ° exp j- ^ dg ^ - ^ } , (12) 
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where for the anomalous dimension function 7(g) and the /3-function the ex- 
pressions up to three loops may be inserted for values of g small enough that 
perturbation theory can be trusted: 

i(g 2 (ri) = iog 2 (v) + W(/i) + 720 V), (13) 



(3(g\ri) = /VM + Pig\v) + /W(aO- (14) 

We note that for 7(g) we know the effective three-loop term from our non- 
perturbative computation of 72 0, while 70 and 71 are given from perturbation 
theory. 

From eq. (§) once the O ren (/^o) is known in some scheme, for example the SF 
scheme we have described, we can obtain the renormalization scheme invariant 
matrix element by introducing an "ultraviolet invariant" running step scaling 
function defined by: 



©inv^o) = *(/////„, f^o)) ■ (9 2 ^)r io/2bo exp { - / dg 



as follows: 



Q ren = Qren^ . 6 UV (/io) . (16) 



i{g) _ 7o_ 

J(g) b g 



(15) 



The scale /io is in general a low-energy scale, where the hadronic matrix element 
can be calculated without severe finite volume effects. In our case, it can be 
identified with a low-energy scale at which the evolution of the renormalization 
constant can be started. In particular we shall fix this scale to be 2L max as in 
ref. [0]. Recently, L max has been computed in terms of the low energy reference 
quantity r in flIU| . In order to "step down" from this scale, we will need the 
step scaling function with s — 2, i.e. starting from ^ 2 (L max ) = 3.48, our largest 
value of g 2 , we evolve with a step size of 2 until contact with perturbation theory 
can be made. 

In this paper, we calculate, as a first step towards the computation of the renor- 
malization group invariant matrix element, the quantity (5|^(/Zo = JIT — )• Note 

2 We remark that the invariance holds with respect to a change of the "ultraviolet" scale /i 
and not of the "infrared" scale /io . 
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that ©J^(yUo) still depends on the reference scale fiQ. The dependence on // will 
only disappear later, when it will be matched with the proper hadronic matrix 
element, making Of^y renomalization scheme independent. 
In order to rely on the perturbative expansion for the (3- and 7-functions appear- 
ing in eq. (|15|), we had to extend the calculation of our non-perturbative running 
to higher scales. We added four more values of g for the step scaling function that 
now covers, in total, values of g 2 (L) ranging from g 2 (L) = 3.48 to g 2 (L) = 0.8873. 
For the results at the four lowest values of g we used the non-perturbatively im- 
proved clover action |TT|; in figs. 1 and 2 we report the continuum extrapolation, 



for the values of g 2 not presented already in ref. || of the step scaling function of 
the quantities fi and fa of eq. (|5|) (a 2 and cr^, respectively, see [[J): at smaller 
length scales the effects of lattice artefacts for 02 are progressively reduced and 
the extrapolations become flatter. 

From the results for o z at the, in total, nine values of g, we can make a fit to the 
step scaling function as a function of g 2 (L). The results for a z at the five largest 
values of g 2 (L) are taken from the combined data presented in ||. In ref. || 
we have shown that, in the scheme we adopted, the coefficient of the two- loop 
anomalous dimension is very large, when compared for example to the one in the 
MS scheme. We have also shown that this coefficient reduces by changing the 
expansion parameter, i.e. by using g 2 (L/4) instead of g 2 {L). The step scaling 
function as a function of g 2 (L/A) is well fitted numerically by a polynomial in 
g 2 (L/4) of the form: 

a(g 2 (L/A)) = 1 - 7 o log(2)^ 2 + c 4 • t + c 6 ■ t + c 8 ■ f , (17) 

where 70 = 4/(97r 2 ). The final results stay unchanged when we switch to a two- 
parameter fit that also gives a very good \ 2 ■ We show our data for o z as a 
function of g 2 (L/4) together with the fit of eq. (|T7|) in fig. 3. 
From this fit we can construct the running step scaling function of eq. (^) with 
fiQ = (2L max )~ 1 . The result is shown in fig. 4, where we have used the two-loop 
expression for 7(g) and the 3-loop expression for (3(g). By using eq. flT5| ) we can 
finally estimate the value of ©^(/io): in the second column of table 1 we report 
the values of ©^(/zo) as a function of the scale [i: for large values of fi the 
function, within the errors, approaches a plateau. We make a fit to a constant 
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Figure 1: Continuum extrapolation of o$ x using a linear fit to the three data 
points with smallest values of a/L for the most perturbative values of g 2 we have 
used in our work, which are indicated in the figure. 
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Figure 2: Continuum extrapolation of 0% using a quadratic fit to all four data 
points for the most perturbative values of g 2 we have used in our work, which are 
indicated in the figure. 
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Table 1: The values for @jj^(/x ) when different scales \i are taken for matching 
with perturbation theory. 

for the results ranging from /j,//j,q = 2 5 to /i//io = 2 9 , and we finally quote: 

6gXr(/io = i^L^)- 1 ) = 1.11(2) . (18) 

The invariant step scaling function is still scheme- dependent, because of the pres- 
ence of the reference scale /xq- This will be cancelled only in the combination that 
defines the invariant matrix element. However, at fixed /Jo, it should be indepen- 
dent of the choice of g 2 (L/A) or of g 2 {L) in the fit to the step scaling function. 
We therefore repeated the whole procedure described above by fitting the step 
scaling function as a function of g 2 (L) and by using the correspondingly modified 
gamma function to two loops. The results are given in the third column of ta- 
ble 1. They are fully compatible with those obtained from the case L/4, although 
the plateau starts at higher energies, as expected. We report the comparison of 
both cases also in fig. ^. In the fourth and fifth column of table 1 we report the 
result for the case "L/4" and "L" respectively, after including our estimate of 
the three-loop anomalous dimensions for the two cases, determined in || U . Not 
surprisingly, the two cases get close to each other more precociously. An estimate 
of the renormalization group invariant yields @J^(/i = — ) = 1-14(2), again 
consistent with our earlier results using g(L/4) as expansion parameter. 
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Figure 3: Our fit to the step scaling function. 
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Figure 4: ©Jj^(/i ) as a function of the scale /x normalized to our reference scale 
/io = (2L max ) 1 . 
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Matching the results of this paper with a non-perturbative calculation of the 
hadronic matrix element, in the continuum, in the same scheme and at the same 
reference energy scale, leads to the definition of a renormalization group invariant 
matrix element that can be confronted with experiment at any scale and in a 
preferred scheme. Such a calculation is in progress. 
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